A test of non-equilibrium thermodynamics in glassy systems: the soft-sphere case 
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The scaling properties of the soft-sphere potential allow the derivation of an exact expression 
for the pressure of a frozen liquid, i.e., the pressure corresponding to configurations which are 
local minima in its multidimensional potential energy landscape. The existence of such a relation 
offers the unique possibility for testing the recently proposed extension of the liquid free energy 
to glassy out-of-equilibrium conditions and the associated expression for the temperature of the 
configurational degrees of freedom. We demonstrate that the non-equilibrium free energy provides 
an exact description of the soft-sphere pressure in glass states. 



The potential energy landscape (PEL) formalism 
has provided a transparent formulation of the equilibrium 
free energy of supercooled liquids based on the statistical 
properties of a system's multidimensional potential en- 
ergy surface |lL|2|,|i2|^yilauLl£|> i- e -j m terms of the 
number, depth, and associated phase-space volume of lo- 
cal potential energy minima . This formalism, in which 
the PEL minima are termed inherent structures (IS), is 
well suited to the description of supercooled liquids and 
glasses. While the liquid explores an exponentially large 
number of distinct PEL basins, a glass explores only a 
very small fraction of these (within an experimentally- 
accessible time frame.) The equilibrium free energy is 
written as sum of two contributions: an entropic term, 
— TS con f(Eis), which accounts for the number of basins 
of depth Eis, and the term E\s + F v n,(T, Bis), which de- 
scribes the free energy of the system confined to an aver- 
age basin of depth Eig. Here, 5 con f is associated with the 
energy degeneracy of mechanically stable configurations 
and hence is termed the configurational entropy, while 
F V ib is associated with the "vibrational," kinetic distor- 
tions of the system around these configurations. Analo- 
gous expressions have been derived from alternative the- 
oretical approaches 0, 0, 0, El . 

The extension of the supercooled liquid free-energy to 
the out-of-equilibrium case has recently been proposed 
[3 Q El El E3, El In the PEL formalism, out-of- 
equilibrium conditions are implemented by imposing the 
constraint that the basin explored by the system while 
aging docs not coincide with the typical equilibrium one 
|l7|. The proposed free energy is given by 

F(V, T, T„) - -T e S COTd (V, E 1S ) + E 1S + F vih (V, T, E 1S ) 

(1) 

The main distinction between this expression and the 
equilibrium case is that while -F v ib is still evaluated at the 



thermostat temperature T, the configurational entropy 
term is weighted by an additional temperature T e , which 
may be thought of as the temperature of the (out-of- 
equilibrium) configurational degrees of freedom. When 
T D and T differ, the system is in a non-equilibrium state. 
While in equilibrium the value of E\s is controlled only 
by T and V via the condition dF(T,V)/dEi S = 0, in 
this non-equilibrium setting the value of Ei$ is also a 
function of T c . In this case the basin depth which the 
system samples, Ejs(V,T,T e ), is the solution of 
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This expression is based on the hypothesis that the out- 
of-equilibrium system samples a distribution of basins 
similar to the one explored in equilibrium |14lll5lll6tll7| . 
Inverting Ejs(V,T, T e ) provides an estimate of T e when 
the bath temperature is T and the system is confined 
to a basin of depth E\s- As already discussed E3> the 
expression for T c coincides with the experimentally de- 
termined Active temperature 

El 

for models in which 
dF v fo/dEis = 0, i.e., when the phase-space volume of 
basins is independent of their depth. The same expres- 
sion for T e has been derived also by Franz and Virasoro 
[16j in the context of p— spin systems, once the basin free 
energy is identified with the TAP free energy |Toj . 

Starting from the proposed free energy (Eq. ^) it is 
possible to calculate the thermodynamic pressure in out- 
of-equilibrium conditions, when the bath temperature is 
T and the configurational temperature is T c (i.e. when 
the system is exploring basins different from those ex- 
plored in equilibrium at temperature T). To this end, 
we evaluate the constant-T volume-derivative of the free 
energy in Eq. ^ 
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and rearranging: 
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The first term in the rhs of Eq. [S] is zero by the Pis 
condition in Eq. [21 Thus, 



P(V,T,T e )=T e 
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The above expression, when evaluated at the bath tem- 
perature T = OK , provides the theoretical expression for 
the pressure experienced in an inherent structure, Pis- 
Indeed, an infinite cooling rate quench to T = OK brings 
the system to the local inherent structure 0,H(J- In other 
words, the steepest descent procedure used numerically 
to locate an inherent structure is equivalent to the phys- 
ical process of setting the bath temperature to T = OK. 
The absolute zero bath temperature eliminates any con- 
tribution arising from the vibrational free energy and 
hence, Pis (V, Eis) is 



Pis(V > E IS ) = P(V,0,T e ) = T e 



Using the mathematical identity 
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and the Pis condition in Eq. [21 evaluated at T = 0, we 
obtain 



Pis (V, Pis) = -T € 
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which provides an alternative definition of the inherent 
structure pressure. 

We now show that the scaling properties of the soft- 
sphere potential provides a consistency check for the de- 
rived expression for Pis and, simultaneously, the valid- 
ity of Eqs. n and [21 Indeed, the thermodynamic ex- 
pression for Pis is explicit in the case of the soft-sphere 



potential, without further assumptions. This potential 
has been extensively studied as a model for liquids and 

glasses n mi m m m m m m m 0. The 

potential energy P of a system composed of N parti- 
cles interacting via a soft-sphere potential is E(r N ) — 



E 



N 



where r 



iV 



with Yi the coordinates of particle 



{ri,r 2 , ....,rjv} 
i, and e and a fix 
the energy and length scales, respectively. 

The self-similar nature of the soft-sphere potential has 
a remarkable property in that uniform scaling of particle 
coordinates does not produce changes the topology of the 
potential energy landscape, since P(Ar N ) = A~ n P(r N ) 
|30L |31| . In the PEL formalism this scaling property im- 
plies that the total number of inherent structures and 
basins is invariant to volume changes. The scaling has 
two important additional consequences which we exploit 
in the present study: (i) an isotropic compression of a 
configuration which is a local potential energy minimum 
remains a local minimum; the potential energy change as- 
sociated with the compression coincides with the change 
in the Pis value; (ii) an isotropic compression moves the 
soft-sphere system along a path of constant configura- 
tional entropy (since the number of basins of depth Pis 
at volume V is identical to the number of basins of depth 
Pis + 8 E\s at volume V + SV). Because Pis is the mea- 
sure of the change of system's energy under a compression 
of an IS*— configuration, these two considerations allow 
us to write Pis as the volume derivative of the inherent 
structure energy along a constant configurational entropy 
path: 



Pis (V, Pis) 
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(10) 



This expression, whose derivation has been based only 
on the self-similar nature of the soft-sphere potential, co- 
incides with the general expression in Eq. [21 derived from 
the out-of-equilibrium thermodynamic approach and pro- 
vides a strong validation of the proposed Eas. 11121 

To the extent that the out-of-equilibrium formalism 
is an appropriate description of systems beyond soft- 
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spheres, the expression for the inherent structure pres- 
sure given by Eq.[5]is quite general. Encouraging support 
for this statement is given by the behavior of the inherent 
structure pressure below an ideal glass transition|32| . 

To summarize, the present Letter provides support 
for the recently proposed out-of-equilibrium approach to 
supercooled liquids and its corresponding definition of 
configurational temperature (i.e., the temperature char- 
acterizing a system's sampling of inherent structures 
which are distinct from those it would sample in equi- 
librium.) In the context of the soft-sphere system, this 
non-equilibrium formalism provides a consistent expres- 
sion for the inherent structure pressure, Pis- Further- 
more, the results presented in this Letter provide a for- 
mal and general derivation of Pis (Eq. GJ) in terms of 
statistical properties of the landscape |g, Il8| and open 
the way for a consistent formulation of thermodynamic 
properties in disordered materials based on a separation 
of configurational and vibrational properties. 
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